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Abstract. An empirical calibration is presented for the synthetic Lick indices (e.g. Mgj, (Fe), H/?, etc.) of Simple Stellar 
Population (SSP) models that for the first time extends up to solar metallicity. This is accomplished by means of a sample of 
Milky Way globular clusters (GCs) whose metallicities range from ~ Zq/30 to Z ~ Zq, thanks to the inclusion of several 
metal rich clusters belonging to the Galactic bulge (e.g., NGC 6553 and NGC 6528). This metallicity range approaches the 
regime that is relevant for the interpretation of the integrated spectra of elliptical galaxies. It is shown that the spectra of both 
the globular clusters and the Galactic bulge follow the same correlation between magnesium and iron indices that extends 
to elliptical galaxies, showing weaker iron indices at given magnesium indices with respect to the predictions of models that 
assume solar-scaled abundances. This similarity provides robust empirical evidence for enhanced [a/Fe] ratios in the stellar 
populations of elliptical galaxies, since the globular clusters are independently known to have enhanced [a/Fe] ratios from 
spectroscopy of individual stars. We check the uniqueness of this a-overabundance solution by exploring the whole range 
of model ingredients and parameters, i.e. fitting functions, stellar tracks, and the initial mass function (IMF). We argue that 
the standard models (meant for solar abundance ratios) succeed in reproducing the Mg-Fe correlation at low metallicities 
([Z/H] < — 0.7) because the stellar templates used in the synthesis are Galactic halo stars that actually are a-enhanced. The 
same models, however, fail to predict the observed Mg-Fe pattern at higher metallicities ([Z/H] ^ —0.7) (i.e., for bulge clusters 
and ellipticals alike) because the high-metallicity templates are disk stars that are not a-enhanced. We show that the new set of 
SSP models which incorporates the dependence on the [a/Fe] ratio (Thomas, Maraston & Bender 2002) is able to reproduce 
the Mg and Fe indices of GCs at all metallicities, with an a-enhancement [a/Fe] = +0.3, in agreement with the available 
spectroscopic determinations. The H/3 index and the higher-order Balmer indices are well calibrated, provided the appropriate 
morphology of the Horizontal Branch is taken into account. In particular, the Balmer line indices of the two metal rich clusters 
NGC 6388 and NGC 6441, which are known to exhibit a tail of warm Horizontal Branch stars, are well reproduced. Finally, we 
note that the Mg indices of very metal-poor ([Z/H] < — 1.8) populations are dominated by the contribution of the lower Main 
Sequence, hence are strongly affected by the present-day mass function of individual globular clusters, which is known to vary 
from cluster to cluster due to dynamical effects. 

Key words. Galaxy: globular clusters: general - galaxies: ellipticals and lenticular, cD - galaxies: abundances - galaxies: 
evolution - galaxies: formation 



1. Introduction 

Galactic spheroids, i.e., elliptical galaxies and the bulges of 
spirals, include a major fraction (perhaps the majority) of the 
stellar mass in the nearby universe (e.g., Fukugita et al. 1998). 
From the uniformity of their fundamental properties at low as 
well as high redshift {z <,!) it has been inferred that the bulk 
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of stars in spheroids must be very old, likely to have formed 
at z > 2 — 3 (see Renzini 1999 for a comprehensive review; 
see also Peebles 2002). Moreover, the space density of pas- 
sively evolving galaxies at 2 1 is consistent with the bulk 
of massive spheroids being akeady in place and 2-3 Gyr old 
at this early epoch (Cimatti et al. 2002a,b). Current renditions 
of hierarchical galaxy formation in CDM dominated universes 
have so far failed to predict these empirical findings, favoring 



2 



C. Maraston et al.: Lick indices of Bulge stellar populations. II. Implications for SSP models and elliptical galaxies 



instead a late formation with major activity even below z ~ 1. 
The culprit is probably in the ways in which star formation and 
feedback processes have been parameterized and implemented 
in the so-called semi-analytic models of galaxy formation and 
evolution (e.g., Kauffmaim & Chariot 1998; Cole et al. 2001; 
Somerville et al. 2001; Menci et al. 2002). 

Given our poor understanding of the star formation and 
feedback processes, the detailed study of the stellar popula- 
tions in nearby galaxies (the fossil records) can provide im- 
portant clues on their early formation phases, complementary 
to the direct observation of very high-z galaxies. Indeed, the 
ages and metallicities of the stellar populations of a galaxy are 
useful constraints to its formation mechanism. However, the 
determination of absolute ages and metallicities of composite 
stellar populations from their integrated spectra is hampered by 
well known degeneracy effects (Faber 1972; O'Connell 1976; 
Renzini 1986; Worthey 1994; Maraston & Thomas 2000), and 
a further compUcation arises when the abundances of major 
chemical elements (iron, magnesium, oxygen, etc.) are consid- 
ered, as traced by narrow-band spectroscopic indices such as 
the so-called Lick indices, Mg2, (Fe), H/3, etc. (Burstein et al. 
1984; Faber et al. 1985; Worthey et al. 1994). The appHcation 
of this technique to ellipticals revealed that the observed rela- 
tion between magnesium and iron indices (Mg2 or Mg b vs. 
(Fe)) disagrees with the predictions of population synthesis 
models where the Mg/Fe ratio is assumed to be solar. Observed 
magnesium indices at given iron are significantly stronger than 
in the models (see Fig. 2 in Worthey et al. 1992). This finding 
was confirmed by several subsequent studies for other samples 
of ellipticals (Davies et al. 1993; CaroUo & Danziger 1994; 
Fisher et al. 1995; Jorgensen 1997, Kuntschner & Davies 1998; 
Mehlert et al. 1998; Longhetti et al. 2000; Thomas et al. 2002a.) 

If Lick indices of magnesium and iron trace the correspond- 
ing element abundances, and the models that are meant for 
solar ratios of these elements are correct, then the observed 
indices imply a supersolar Mg/Fe ratio in ellipticals. In turn, 
according to common wisdom this implies short {t<^l Gyr) 
star formation timescales for the stellar populations of ellip- 
ticals (e.g., Matteucci 1994; Thomas et al. 1999). In fact, the 
so-called a-elements (i.e., O, Mg, Ca, Ti, and Si) are promptly 
released by massive, short-living ( ^ 3 x 10^ yrs) progeni- 
tors exploding as Type 11 supernovae, while most iron comes 
from Type la supernovae, whose progenitors span evolutionary 
timescales from over ^ 3 x 10*^ yrs to many Gyrs (e.g., Greggio 
& Renzini 1983; Matteucci & Greggio 1986; Pagel 2001)). 
Therefore, a high a-over-iron ratio ([a/Fe]) implies that star 
formation ceased before the bulk of Type la supernovae had 
the time to enrich with iron the interstellar medium while this 
was still actively forming stars. Such a short star formation 
timescale appears to be at variance with the predictions of cur- 
rent hierarchical models for the formation of elliptical galax- 
ies (Thomas 1999; Thomas & Kauffmann 1999), which predict 
star formation to continue for several Gyrs. In conclusion, the 
Lick indices of magnesium and iron appear to offer a unique 
opportunity to estimate the timescale of star formation in galac- 
tic spheroids, hence to help for a better understanding of the 
early formation phases of galaxies. 



However, a caveat is in order over the above chain of argu- 
ments: how well do Lick indices trace element abundances? 
Are we sure that population synthesis models correctly pre- 
dict the values of these indices? (for early discussions of these 
issues see Tripicco & Bell 1995; Greggio 1997; Tantalo et 
al. 1998). Indeed, the population synthesis models on the ba- 
sis of which the magnesium overabundance has been inferred 
were not calibrated, especially in the metallicity range (so- 
lar and above) which is relevant to elliptical galaxies. Hence, 
it could not be excluded that the population synthesis mod- 
els would underpredict the strength of the magnesium indices, 
while the Mg/Fe ratio of ellipticals would actually be solar By 
calibration of the indices we mean the comparison of their syn- 
thetic values with the corresponding quantities measured on 
objects for which the age and the detailed chemical compo- 
sition - total metallicity and element abundance ratios - are in- 
dependently known. For this comparison the best stellar pop- 
ulation templates are the Galactic globular clusters. However, 
existing databases of Lick indices of globular clusters (Burstein 
et al. 1984; Covino et al. 1995) are restricted to the metal-poor 
objects of the Halo. The most metal rich cluster in the Covino 

et al. sample is 47 Tuc ([Fe/H] 0.7) whose Mg^ is ~ 0.18 

mag, much less than found among elUpticals, which span from 
Mg2~0.2 to-- 0.4 mag. 

Globular clusters that are more metal rich than 47 Tuc 
do actually exist in the Galactic bulge, reaching Z ^ Zq 
for NGC 6553 and NGC 6528 (Barbuy et al. 1999; Cohen 
et al. 1999). Cohen et al. (1998) have measured some of 
the Lick indices for these two clusters, but did so using the 
Burstein et al. (1984) passbands to define the indices. This 
does not allow a direct comparison with the theoretical mod- 
els, which are based on the passbands defined by Worthey et 
al. (1994). Gregg (1994) measures and analyses spectral in- 
dices for several Milky Way GCs including metal-rich objects, 
among which NGC 6528. Though similar, these spectral in- 
dices are not in the Lick system. In Sect. 3 we show that in 
a model calibration it is crucial that data and models are set 
up on the same system, because there are sizable differences 
in the value of some of the indices, depending on the adopted 
passbands. Therefore, we obtained optical spectra for a sample 
of Bulge globular clusters with metallicities [Fe/H] ^ — 0.5 
(including NGC 6528 and NGC 6553), plus some metal-poor 
globular clusters in order to check the models on a wide metal- 
licity range. In fact, existing Lick indices of metal poor globular 
clusters (e.g. Covino et al. 1995: Cohen et al. 1998) were also 
measured in the Burstein et al. (1984) system. The results of 
the measurement of the indices in the Lick/IDS system are re- 
ported in an accompanying paper (Puzia et al. 2002, hereafter 
Paper I). 

For at least some of the program clusters the abundance 
ratios of a-elements to iron are known from high resolution 
spectroscopy of individual stars in these clusters (Barbuy et 
al. 1999; Cohen et al. 1999, Carretta et al. 2001; Coelho et 
al. 2001). While there is certainly room for further improve- 
ments in the abundance determinations, these studies indi- 
cate an overabundance [a/Fe] ~ 0.2 — 0.3 for these clus- 
ters. Moreover, their age, determined from color-magnitude 
diagrams, is virtually identical to the age of Halo globular 
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clusters, i.e., 12-13 Gyr (Ortolani et al. 1995; Rosenberg et 
al. 1999; Feltzing & Gilmore 2000, Zoccali et al. 2001). Having 
fairly accurate estimates for their basic parameters (age, [Fe/H], 
[a/Fe], the measured Lick indices of Bulge globular clusters 
are used to calibrate the population synthesis models and to 
test the "magnesium overabundance" solution for ellipticals. 

The paper is organized as follows. A brief summary of the 
database is presented in Sect. 2. In Sect. 3 a critical overview of 
existing ambiguities in the definition of model and data metal- 
licities, is given. The main results are presented in Sect. 4, 
where the Lick indices of the Halo and Bulge clusters are com- 
pared to those of elliptical galaxies. The following Sections 
explore in detail several technical aspects of the population 
synthesis modeling (Sect. 5) and the calibration of the models 
(Sect. 6). Readers mainly interested in the implications of the 
results on galaxies can skip these sections. Sect. 7 summarizes 
the conclusions. 



2. The Globular Cluster Data 

Optical spectra (3400 < A < 7500 A) have been obtained with 
the Boiler & Chivens spectrograph at the ESO 1.5m telescope 
for a sample of 12 globular clusters (GCs) mostly located in 
the Galactic bulge and for several (15) positions in the bulge 
field known as Baade's Window. The data acquisition, reduc- 
tion and the resulting indices in the Lick system are fully de- 
scribed in Paper I. Here we summarize some key features of 
the data, which are useful to the present discussion. 

The target GCs have been selected on the basis of two 
requirements. First, a high metallicity in order to extend the 
model calibration towards the range most relevant to ellipticals. 
Among the 12 clusters in the sample, 7 clusters have metallic- 
ities [Fe/H] ^ - 0.5 (on the scale of Zinn & West, 1984), the 
most metal-rich ones being NGC 6528 and NGC 6553. The 
remaining 5 clusters are more metal-poor, and were included 
in the sample to check consistency with previous studies (e.g. 
Trager et al. 1998). Second, the availability of independent es- 
timates of element abundances, total metaUicities and ages, in 
order to allow for the empirical calibration of the synthetic in- 
dices. As already mentioned, estimates of the metallicity and 
[a/Fe] ratios are available for the two well-studied clusters 
NGC 6553 and NGC 6528, which ensures a meaningful model 
calibration around solar metallicity. For the remaining clusters 
estimates of the metalhcity in the Ziim & West scale and in the 
Carretta & Gratton (1997) scale are available. This allows us to 
use the clusters to calibrate the metallicity scale of the models 
in a relative sense. Nevertheless, it would clearly be useful to 
extend the detailed elemental abundance determinations to all 
the clusters in the sample. 

Special care has been paid to subtract the foreground/back- 
ground Ught from the cluster's Ught. Indeed, the field and 
clusters stellar population components in the Bulge appear 
to be virtually coeval, and to span similar metallicity ranges 
(Ortolani et al. 1995; Zoccali et al. 2002). The very similar, 
though not identical, spectral energy distribution of the bulge 
light could therefore introduce spurious effects on the measure- 
ment of the cluster indices, if not adequately subtracted. 



The luminosities sampled in the GCs and the bulge fields 
by the slit of the spectrograph have been carefully evaluated 
in order to assess the dependence of the indices on stochastic 
effects. The number of stars that are expected to be detected 
in the various evolutionary phases is proportional to the total 
sampled luminosity and can be easily evaluated (Renzini & 
Buzzoni 1986, Maraston 1998, Renzini 1998). Therefore, for 
every cluster it has been checked whether the sampled lumi- 
nosity is dominated by few, very bright stars, like RGB-tip or 
E-AGB stars (Table 2, Paper 1). This is important for metal- 
lic Lick indices like Mgj, TiO, NaD, which are very strong in 
these stars. The uncertainties on the indices associated with the 
stocastic fluctuation in the number of stars which contribute the 
light in the relevant wavelength ranges, are included in the error 
budget (see Paper 1). 

3. A note on the metallicity calibration of SSP 
models 

Some of the Lick indices were designed as metallicity indica- 
tors for unresolved stellar populations, therefore to calibrate a 
model Lick index means to check whether the model gives the 
observed value of the index for a SSP whose age and compo- 
sition are independently known. In practice, GCs offer the best 
proxy to a SSP. However, some ambiguities make such calibra- 
tion not so straightforward. 

3. 1. Ambiguities in the definition of metallicity 

From the model side, the total metallicity of model Lick in- 
dices is not well defined because of the role of the so-called 
"fitting functions" (see Sect. 5.1). The fitting functions are 
best fits of the Lick indices as measured in stars, as func- 
tions of the stellar parameters Tcs, g and chemical compo- 
sition. According to the standard procedure (e.g. Buzzoni et 
al. 1992, 1994; Worthey 1994), the fitting functions are plugged 
on the isochrones to compute the Lick indices of SSP mod- 
els (Sect. 5, Eqs. 4-5). Therefore, it is necessary to specify 
the metallicity parameter(s) for the fitting functions and for the 
isochrones, and of course they should be the same. While the 
latter is well-defined by construction (stellar evolutionary mod- 
els are constructed for well defined sets of abundances), the 
former is somewhat ambiguous. Indeed, the estimates of the 
chemical composition of the stars used for the fitting functions 
come from a variety of sources, both spectroscopic and photo- 
metric. These, quite inhomogeneous, metaUicities are collected 
under a parameter referred to as [Fe/H] in the fitting func- 
tions available in the Uterature (Worthey et al. 1994; Buzzoni 
et al. 1992; 1994). An additional source of compUcation comes 
from the fact that a certain total metallicity might be achieved 
with different proportions of the major elements, the so-called 
a-elements, with respect to iron. The fact that the fitting func- 
tions are derived from observed stars implies that the element 
abundance ratios are not constant in the fitting. In fact, as 
well known, the [a/Fe] ratios vary systematically among Milky 
Way stars (e.g. Mc William 1997), including those in the sam- 
ples used to construct the fitting functions themselves. On the 
other hand the specific abundances of magnesium and iron, be- 
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Table 1. Lick indices of 47 Tuc (spectra of Covino et 
aL 1995) measured in the Worthey et al 1994 (W94) sys- 
tem for this work. The corresponding values in the Burstein 
et aL 1984 (B84) system by Covino et aL (1995) are given in 
second line. 



Lick System 


Mg2 


Mg& 


Fe5270 


Fe5335 


H/3 




(mag) 


(A) 


(A) 


(A) 


(A) 


W94 


0.15 


2.19 


1.64 


1.40 


1.48 


B84 


0.18 


3.02 


2.18 


1.88 


1.62 



side the total metallicity likely affect the strength of Mg and Fe 
absorption lines (see Sect. 4). 

From the GC data side, the empirical metallicity scale of 
GCs is not rigorously defined either. The reference values for 
the chemical composition of the sample GCs used in this work 
are taken from the revised compilation by Harris (1996), which 
is largely based on the Zinn & West (1984) scale. Thomas et 
al. (2002b) show that the Zinn & West (1984) scale, which is 
named as [Fe/H], is likely to be closer to the total metallicity 
rather than to the sole iron abundance. In fact, the Zinn & West 
scale is tied to the scale set up by Cohen (1983), in which the 
metallicities which are called [Fe/H] are indeed obtained by 
averaging [Mg/H] and [Fe/H] (Thomas et al. 2002b). This fact 
was anticipated by the evidence that the integrated colours of 
SSP models as function of the model total metallicity match 
well with those of Milky Way GCs, when the metallicities of 
the latter are on the Zinn & West scale (Maraston 2000, Fig. 1). 
Moreover, the values reported in the Harris (1996) catalogue 
are «of just the metallicity in the Zinn & West scale in all cases. 
When various estimates of metallicity from other sources are 
available, either from spectroscopy or from colour magnitude 
diagrams, these are used together with the Zinn & West-based 
value, and the straight average of the values is published as 
[Fe/H]. 

3.2. The effect of the adopted Lick system 

A quantitative model calibration requires that data and mod- 
els refer to the same spectro-photometric system (see Maraston 
et al. 2001a). The definition of the Lick system (index pass- 
bands, resolution, etc.) has been slightly changed from Burstein 
et al. (1984), through Worthey et al. (1994, hereafter W94) to 
Trager et al. (1998). The different index definitions may in- 
troduce offsets, which could affect the model calibration. We 
adopt here the data as measured in the W94 version of the Lick 
system, since the models are locked to this version via the in- 
dex fitting functions. The effect of this choice has been tested 
by computing the W94-like indices for the GC spectra obtained 
by Covino et al. (1995), then comparing them to the Covino 
et al. values of the indices that are measured in the Burstein 
et al. (1984) system. The comparison for the case of 47 Tuc 
is shown in Table 1. The differences among the indices are 
not negligible. The values on the W94 system are systemati- 
cally lower than those on the Burstein et al. (84). In Sect. |^ we 
demonstrate the impact of the index definitions on the model 
calibration. 



4. Results 

Fig. 1 shows the Mg b vs. the average iron index (Fe) [] 
of the GCs of our sample (large filled symbols). The large 
open circle refers to the coaddition of the spectra of 15 Bulge 
fields located in Baade's Window. The blue lines are models 
of Simple Stellar Populations (SSPs), i.e. coeval and chemi- 
cally homogeneous stellar populations, with total metallicities 
([Z/H]) Q ranging from —2.25 to +0.67, and ages between 3 
and 15 Gyr. The models are computed with the evolutionary 
population synthesis code of Maraston (1998), as described in 
Maraston & Thomas (2000; see also Maraston et al. 2001b 
and Maraston 2002), and are based on stellar tracks, imple- 
mented with the Worthey et al. (1994) fitting functions, to de- 
scribe the stellar indices as functions of effective temperature, 
gravity and metallicity. The stellar evolutionary tracks are from 
Bono et al. (1997) and Cassisi et al. (1999) for metallicities up 
to Z = 0.04, from Salasnich et al. (2000) for Z = 0.07, and 
adopt solar abundance ratios. In the following, we refer to these 
models as standard SSPs. ^ 

The GCs indices define a nice sequence to which the Bulge 
field appears to belong as well. The sequence runs with a shal- 
lower slope compared to the standard models, i.e. at a given 
(Fe) index the data have a stronger Mg 6 than the models. 
Several observational evidences show that all galactic GCs 
have a enhanced abundance ratios, with typical values around 
[a/Ye] ~ +0.3 (e.g. Pilachowski et al. 1983; Gratton, 1987; 
Gratton & Ortolani 1989, Carney 1996, Salaris & Cassisi 
1996). In particular, for the two most metal rich clusters (NGC 
6553 and NGC 6228 with Mg - 3.8) Barbuy et al. (1999) find 
[a/Fe] ^ +0.3 from individual star spectroscopy. The Bulge 
field stars are also known to be overabundant in Mg with re- 
spect to the solar ratio (McWilliam and Rich 1994). 

Therefore the GC sequence traces the locus of a-enhanced 
SSPs. 

Fig. |l]aIso shows the central values (i.e. those obtained with 
apertures R ^ 1/8 + 1/10 Re) of the indices of field and 
cluster ellipticals taken from various sources in the literature 
(see caption). The indices of ellipticals occupy a relatively nar- 
row range in (Fe) and a large range in Mg 6, stretching from 
the standard models to the high metallicity extension of the 
GCs sequence, and beyond. With very few exceptions, both 
Mg h and (Fe) indices are measured stronger in the nuclei of 

' (Fe)=(Fe5270+Fe5335)/2 

" The notation [Z/H] is used to indicate total metallicities, i.e. the 
total abundance of heavy elements with respect to hydrogen nor- 
malized to the solar values, i.e. [Z/H] — log(-^) — log( jS^). By 
[Fe/H] we mean the abundance of iron with respect to hydrogen nor- 
malized to the solar values, i.e. [Fe/H] — log(T7^) — log(f|3-). If 

elements have solar proportions then [Fe/H] = [Z/H]. In case of 
«-element enhancement, the relation between [Fe/H] and [Z/H] is: 
[Fe/H] = [Z/H] - 0.94 * ([a/Fe]) (Thomas et al. 2002b; see also 
Trager et al. 2000). 

' We want to emphasise that what we call standard SSPs, i.e. those 
based on the Worthey et al. or on the Buzzoni et al. fitting functions are 
not solar-scaled SSPs at every metallicity. Indeed this type of models 
are constructed by adopting the stellar indices of Milky Way stars, 
which have a variety of abundance ratios, see Sect. 4. 
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Fig. 1. Mg b indices vs. average iron (Fe) indices of galactic and Bulge GCs of our sample (large red filled symbols). The open 
circle shows the average value of 15 Bulge fields located in Baade's window. Small open grey symbols show the central values 
of the indices for ellipticals taken from the literature: the field and cluster ellipticals from Beuing et al. (2002, pentagons); the 
Coma ellipticals from Mehlert et al. (2000, stars); the Fornax ellipticals of Kuntschner & Davies (1998, squares); the field and 
Virgo ellipticals of Gonzalez (1993, triangles). Standai'd SSP models with metallicities [Z/H] = (-2.25, - 1.35, - 0.55, - 
0.33; 0.00; + 0.35; +0.67), from bottom to top, and ages 3, 5, 10, 12 and 15 Gyr, from left to right, are shown as blue lines. 
The thick black lines show 12-Gyr SSP models with same total metallicities, and various [a/Fe] = 0, + 0.3, + 0.5 (Thomas et 
al. 2002b). 



ellipticals than in the most metal rich GCs in our sample. The 
stellar populations in the centers of ellipticals seem to be char- 
acterized by: 

(i) a supersolar total metalUcity; 

(ii) a range in abundance ratios, from almost solar, to 
[a/Fe] values as large as those of the most metal rich bulge 
clusters, or even more. 

Similar conclusions have been proposed in the literature 
(Worthey et al. 1992; references in Introduction). Yet, these 
were based on the assumptions that standard models reproduce 
the indices of solar abundance ratios SSPs (at least at solar 
metallicity and above), that Mg b and (Fe) trace the Mg and 
Fe abundance, and that an a- element overabundance affects 
the indices in the appropriate direction. 

Our comparison of the GCs data with the ellipticals, and 
especially the inclusion of the indices of the metal rich clusters 



in the galactic Bulge and the bulge field, confirm the validity of 
these assumptions, the clusters being used as empirical SSPs. 

This empirical evidence motivates the construction of new 
SSP models with various, well-defined [a/Fe] ratios which are 
shown in Fig. |l] as thick black lines (Thomas et al. 2002b, here- 
after TMB02). The GCs are now very well represented by a co- 
eval (12 Gyr old) sequence of models with various metallicities 
and [a/Fe] — +0.3, in agreement with the results from stel- 
lar spectroscopy. Note that at low metallicities (Mg 6 < 2) the 
standard models (blue lines) match with the enhanced models. 
This is due to the standard calibrations by Worthey et al. (1994) 
being a-enhanced at sub-solar metallicities (Sect. 5; see also 
TMB02). 

In the next Sect, we conduct a thorough analysis of the stan- 
dard SSP models, with the aim of assessing whether effects 
other than an enhanced [a/Fe] ratio can explain the deviation 
of the data from the standard models. In other words, we in- 
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vestigate the uniqueness of the "magnesium overabundance" 
solution. 

5. Model Lick indices: Icey ingredients and 
ambiguities 

In the model SSPs, the line-strength of an absorption line with 
bandpass A /^sp jg given by 

/ = A-(l-^^i/^^c), (1) 

where F\ and Fc are the fluxes in the absorption line and in the 
continuum, respectively. 

In case of Lick indices this formula cannot be applied 
straightforwardly because of the different spectral resolu- 
tion of the Lick system (~ 8 A) and of the model atmo- 
spheres (Kurucz-based stellar atmospheres have resolutions of 
20 A in the optical region where the Lick indices are de- 
fined). Basically, the problematic quantity is the flux in the 
absorption line F\ , because it depends on the spectral res- 
olution. To overcome this problem, the Lick group has mea- 
sured the Lick indices on observed stellar spectra having the 
required resolution (Burstein et al. 1984; Faber et al. 1985). 
Assigning to each star of their sample the values for the stel- 
lar parameters surface gravity (g), effective temperature (Tcs) 
and chemical composition, they have constructed polynomial 
best-fitting functions which describe the various Lick indices 
measured on the stars, /*, as a function of these parameters, 
i.e. /* — i{Tctf;g; [Fe/H]). These polynomial fittings follow- 
ing Gorgas et al. (1993) are called fitting functions (hereafter 
FFs). 

The integrated Lick index of an SSP model is then evalu- 
ated as it follows. 

The flux in the absorption line of the generic i-th star of the 
SSP, F*j , can be expressed with Eq. 1 as: 

Fl ^ Ki ■ (1 - (2) 

where /* that is the index of the i-th star is computed by in- 
serting in the FFs the values of {Tcs;g; chemical composition) 
of the i-th star, and F* ; is its continuum flux. The latter is 
computed by linearly interpolating to the central wavelength 
of the absorption line, the fluxes at the midpoints of the red and 
blue pseudocontinua flanked to the line (Table 1 in Worthey et 
al. 1994). Eq. 1 can be re-written as: 

/SSP ^ A . (1 - ^ FyJ2 Kd = E ■ fc* i (3) 

i i i 

where /*j are the contribution of each individual star to the 
total continuum flux of the SSP. Thus, the SSP integrated in- 
dex /SSP is the weighted average of the stellar indices /* with 
the weigths being /*;. When computing actual models, the 
isochrone representing the SSP is binned in Tcff subphases, 
small enough to ensure that /j* is the same for the stars be- 
longing to the given subphase. A good binning is ATcff ~ 
100 K (Maraston 1998). Eq. 3 can re-expressed for the sub- 
phases j 

i^'^-T^^r-f:, (4) 



where /*j = Eie, /cV 

It should be noted that since the role of the stellar continua 

is that of a weight, it is not crucial that they are evaluated on 
Kurucz-type spectra. A relation similar to Eq. 4 holds for in- 
dices measured in magnitudes (e.g. Mg2): 

10-o.4.Mgf - ^ J2 io-o-4 Mg^ . J* , 
j 

The two ingredients (/j* and f*j) are discussed comprehen- 
sively in the following subsections. 

5. 1. Interplay between fitting functions rand 
continua 

To explore the systematic effects introduced in SSP models 
by the use of different sets of FFs, and following Maraston 
et al. (2001b), we compute the same SSP models with three 
formulations for the FFs from the literature, i.e. by Worthey et 
al. (1994, hereafter Worthey et al. FFs), Buzzoni et al. (1992, 
1994, hereafter Buzzoni et al. FFs), and Borges et al. (1995, 
hereafter Borges et al. FFs). 

Worthey et al. FFs^ the most widely used in the SSP mod- 
els in the literature, are based on the Lick sample of ~ 400 
nearby stars. Buzzoni et al. FFs are based on a smaller sam- 
ple of stars 87), also located in the solar vicinity. As is well 
known, the a to Fe abundance ratios in nearby stars vary with 
metallicity, ranging from the super-solar values in the Halo 
stars ([ck/Fe] ^ +0.3) to the solar proportions in the metal-rich 
disk stars ([a/Fe] = 0.0) (e.g. Wheeler et al. 1989; Edwardsson 
et al. 1993; Fuhrmann et al. 1995; Fuhrmann, 1998; see the 
comprehensive review by Mc William 1997). Thus, likely, these 
two sets of FFs reflect a enhanced mixtures at low Z, and solar 
ratios at high metallicity, a trend which is dragged into the SSP 
models through /j* (Eqs. 4-5). This explains why the standard 
models in Fig. |l] (blue lines) represent well the indices of metal 
poor GCs. 

Borges et al. FFs (see also Idiart & de Freitas 
Pacheco 1995) have been derived from a sample of roughly 
90 stars for which the Mg to Fe ratio has been measured. Thus, 
they include the [Mg/Fe] ratio as an additional variable besides 
temperature, gravity and metallicity. 

In the following subsections we discuss separately high 
metallicity (Sect. 5.1.1) and low metallicity (Sect. 5.1.2) model 
indices. 

5.1.1. The metal-rich zone 

Fig. ^illustrates the interplay between /j * and /* j (Equation 4- 
5) in determining the SSP magnesium and iron indices in the 
particular case of a 15 Gyr, solar metallicity and solar abun- 
dance ratio SSP, with Salpeter IMF. The three sets of FFs con- 
sidered are color coded as marked in the top-right panel. The 
X-axis is a monothonic coordinate along the SSP isochrone. 
The integer values of a; ( 1 to 6) mark the end of the six main 

We notice a typo in Table 3 of Worthey et al. 1994. The fitting 
functions are mistakenly given as function of log 9, while they are 
function of 6, as also stated in the text. 
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Fig. 2. Closer look to model Mg2 (left-hand panels) and (Fe) (right-hand panels) indices of a 15 Gyr old SSP with solar metal- 
licity, solar abundance ratios, and Salpeter IMF. The upper panels show the cumulative SSP indices /^sp (^^^^ integrated up to the 
point X along the isochrone (see text): x = 1: lower Main Sequence (d-MS) up to Tes < 5000 K; x ~ 2: Main Sequence up to 
Turn Off (TO); x = 3: Sub Giant Branch (SGB); x = 4: Red Giant Branch (RGB); x = 5: Horizontal Branch (HB); x = 6: Early 
Asymptotic Giant Branch (E-AGB). These phases are identified on the corresponding isochrone (small inserted panel), shown in 
the temperature-gravity plane. Here the label [a/Fe] — refers to the composition of the stellar tracks. Colors code the different 
FFs for /j* (Eq. 4-5) adopted for the SSP models, by Worthey et al. (red lines), Buzzoni et al. (green lines) and Borges et al. 
(blue lines). The values /j* are shown in the middle panels, with lines connecting the stellar indices of the subphases along the 
isochrone. Finally the lower panels display the contributions of the various subphases to the total continuum flux, /* j (Eq. 4-5). 
Note that the flux contributions in post-MS (x > 2) take into account the fuel consumption, i.e. the product of the stellar lifetimes 
with the stellar luminosities. For example the contribution of the long lasting RGB bump phase (peak in phase 3-4) is much 
larger than that of the RGB-tip (end of phase 4), where, in spite of larger stellar luminosities, the evolutionary timescale is much 
shorter. 
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evolutionary phases (see the caption and the figure with the 
isochrone inserted in the top-left panel). Each x-point in Main 
Sequence (x < 2) represents the subphase along the Main 
Sequence isochrone. Each z-point in post-MS (x > 2) rep- 
resents the j subphases of every post-MS major phase, which 
are equally spaced in effective temperature (ATcff ~ 100 K , 
Maraston 1998). The top panels show the cumulative Mg2(x)[] 
(left) and (Fe)(a;) (right) along the isochrone, which assumes 
the value of the SSP model at a; = 6. As in Eq. 4-5, each value 
of Mg2(x) and (Fe)(2:) is obtained by summing up to point 
X, the product of /j * times /* j . These are separately shown 
in the central and lower panels, respectively. The behavior of 
/j * in the central panels reflects the changing of Toff along the 
isochrone, both Mg2 and(Fe) indices being very strong in cool 
stars. 

The three sets of EEs correspond to quite different val- 
ues for /* along the isochrone, particularly in the faint dwarf 
regime (x <i 1), at the Tip of the RGB {x ~ 4), and at the end 
of the E-AGB (x ^ 6). In spite of that, the SSP indices keep 
very close along the isochrone, and assume quite similar total 
values, due to the low contribution to the total continuum flux 
of these particular phases. Thus, the indices turn out to be quite 
insensitive to the adopted set of EEs. 

As a result of the weighting through the lower MS, 
RGB and E-AGB bright stars (which have very strong Mgj) 
are not important in determining the total SSP index. The most 
important contributors to the continuum fluxes in the two win- 
dows of Mg2 and (Ee) are: the stars around the TO; those on the 
fainter portion of the RGB (especially at the so-called bump), 
and the Horizontal Branch stars. Their indices dominate the in- 
tegrated values. This applies in general to old ( ^ 3 Gyr) stellar 
populations, as the relative contribution to the total optical flux 
of the different phases does not depend much on age in this age 
range (Renzini & Buzzoni 1986, Maraston 1998). 

The contributions to the total continua of Mg2 and(Ee) of 
the various evolutionary phases are given in Table 2 for 15 Gyr 
old metal-rich and metal-poor SSPs. 

As apparent in Eig. |[ Borges et al. EEs provide /j* for 
faint dwarfs which are much larger than those of the other two 
sets. These result from the exponential increase with decreas- 
ing Tcff of their EE for Mg2. In the validity range as specified 
by the authors (Tes ^ 3800 K), the EE for Mg2 yields values 
as high as 1.9 mag (corresponding to 0.83 A in the units of 
Eig. while the coolest dwarf in their observed sample has 
Mg2 = 0.45 mag (or 0.34 A). These very strong (extrapolated) 
Mg2 indices as obtained with a blind use of the EEs are ex- 
tremely unrealistic. This example illustrates the importance of 
checking the behaviour of the algebraic EEs when computing 
SSP models. 

We conclude that the Mg-Ee relation of standard models 
around solar metallicity (Eig. [I]) is independent of the fitting 
functions. 



^ Notice that we plot the Mgj index expressed in A, i.e. Mgj (A) = 
1 - 10(-'' '' '^g2 (mag)) ]y[g_^ instead of Mg b because 

neither Buzzoni, nor Borges give FFs for Mg b. The two indices, how- 
ever, are very closely related. 



Table 2. Relative contributions of stellar evo- 
lutionary phases, to the continuum flux of the 
Mg2 index (A - 5175A) and (Ee) (A ~ 
5300A), for 15-Gyr old SSPs with Salpeter 
IME, and metaflicities: Zq and Z = lO""'. 
RGB tip is the portion within 1 mag from the 
tip. 



Zq IQ-^Zq 



phase 


Mg2 


(Fe> 


Mg2 


(Fe> 


d-MS 


0.11 


0.13 


0.10 


0.10 


MS TO 


0.32 


0.22 


0.26 


0.26 


SGB 


0.05 


0.12 


0.10 


0.10 


RGB 


0.29 


0.30 


0.28 


0.28 


RGB-tip 


0.01 


0.01 


0.10 


0.10 


HB 


0.17 


0.17 


0.12 


0.12 


E-AGB 


0.05 


0.05 


0.04 


0.04 



5.1 .2. The metal-poor zone 

As shown in the previous paragraph, the differences in the EEs 
appear to be unimportant at ^ Zq. Actually, the Mg2 index 
obtained with the Borges et al. EEs does deviate from the other 
two values, by an amount which is comparable to the typical 
observational error affecting the GC data 0.01). This dis- 
crepancy becomes more pronounced at very low metallicities, 
which is relevant for the model calibration with GCs. 

Eig. 1^ is the analogous of Eig. ^, but for a very metal-poor 
SSP with metallicity Z = 10^^. We use here the Borges et al. 
EEs with [Mg/Ee]=0.3 which is appropriate at low metallicities. 

At low metallicity the temperature distribution of the iso- 
chrone shifts to hotter values (see inserted panels in Eig.s ^ 
and H). The relative contributions to the continuum flux of the 
different evolutionary phases is very similar to the Zq case 
(see Table 2), but the stellar indices are now very weak, and 
the strongest Mg2 are found on the lower MS, which is the 
coolest portion of the isochrone. As a result, the total Mg2 of 
the SSP is very close to the value attained already at the turn- 
off point. Borges et al. EEs yield much lower stellar indices for 
the lower MS, compared to the other two EEs, which explains 
the lower integrated index. Incidentally, this is true also when 
adopting [Mg/Ee]=0. in the EE formula. It should be noted that 
the behaviour of the EEs at the low Main Sequence (phase 1) is 
not constrained by stellar data. Indeed, only 5 main sequence 
stars (symbols in the lower left panel of Fig. ||) are found at 
[Fe/H] w —1.8 by merging both the Lick and the Borges et 
al. data base. An improvement of present low-metallicity SSP 
models can be gained by implementing the stellar libraries with 
cool, dwarf and low metallicity objects. 

The importance of the (lower) main sequence phase on the 
integrated Mg2 has the consequence of making this index sen- 
sitive to the mass function of the stellar population. This com- 
plicates the comparison of the models with the GC data, for 
which the present day mass function derives from the IMF plus 
the possible dynamical evolution, which can lead to the evapo- 
ration of the low mass stars (e.g. Piotto & Zoccali 1999). 
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Fig. 3. Same as Fig. 2 but for a very metal-poor SSP with metal- 
licity Z = 10^^. Note that the integrated Mgj is made up by 
the lower Main Sequence (d-MS, phase 1). The symbols in the 
left-hand bottom panel are the values of Mg2 of the stars in the 
Lick sample (filled) and in the sample used by Borges (open) 
with metallicities and gravities appropriate to the MS of these 
SSP models. As in Fig. ^, the label [a/Fe] = in the inserted 
panel refers to the composition of the stellar tracks. 




0.1 0.2 0.3 0.4 0.5 



Fig. 4. Effect of a dwarf-dominated IMF on the 
Mg2 and(Fe) diagram of SSPs. The model grid for the 
Salpeter IMF exponent of 2.35 is the same as in Fig. |l| The 
two additional solid lines are dwarf-dominated SSPs with 
IMF's exponents of 4 and 5, respectively. These SSPs have 12 
Gyr and metallicities as indicated by the labels. Data of GCs 
and elhpticals as in Fig. |l|. 



Different from the Mgj index, an important contribution 
to the total (Fe) index comes from the RGB portion of the 
isochrone, where the FFs are very noisy (see right panels in 
H). Notice that the index computed with Worthey et al. FFs 
converges to that based on Buzzoni FFs because of negative 
stellar indices predicted for the hot HB stars. It is very difficult 
to assess the reliability of the indices as metallicity indicators 
at such low-Z. 

5.2. IMF effects 

As discussed in the previous sections, Mg2 indices are very 
strong in dwarf stars, therefore dwarf-dominated stellar popu- 
lation could in principle allow to reach the very high values of 
Mg2 shown by ellipticals. Fig. Q shows the location of such 
dwarf-dominated SSPs in the Mg2 vs. (Fe) diagram (lower 
solid lines). Lines connect 12 Gyr old models with metallici- 
ties from 3 time to half solar, and IMF's exponents of 4 and 
5 (in the notation in which Salpeter is 2.35), from top to bot- 
tom, respectively. Worthey et al. FFs have been used for this 
exercise. Data of ellipticals and GCs are the same as in Fig. [l|. 

Dwarf-dominated stellar populations are able to reproduce 
the Mg2 and (Fe) indices of ellipticals without invoking abun- 
dance effects. It should be noted that in such dwarf-dominated 
SSPs, > 70% of the total luminosity is made up by stars close 
to the H-burning hmit. The different slope of these SSP models 



on the Mg2— (Fe) plane reflects the different dependence of 
the stellar indices from Tcff at the lower end of the MS. 

These extreme IMFs are also able to reproduce the low 
values of the Calcium triplet absorption line at 8600 K ob- 
served in ellipticals (Saglia et al. 2002), because its strength 
decreases with increasing gravity (Jones et al. 1984). However, 
these models fail at explaining other spectral properties of el- 
lipticals. In fact the corresponding stellar mass-to-light ratios 
become much larger (M/Lb > 30, see Maraston 1998) 
than the dynamical ones observed in the central portions of 
ellipticals (M/Lb ^ 6, Gerhard et al. 2001). Similarly, a 
dwarf-dominated elliptical galaxy light was excluded given the 
strength of the CO absorption (Frogel et al. 1978) and the ab- 
sence of the Wing Ford bands in absorption (Whitford 1977). 
The case of such a dwarf-dominated present mass function 
for GCs, and the Bulge field is ruled out by direct observa- 
tions of the lower MS (de Marchi & Paresce 1997; Piotto & 
ZoccaH 1999; ZoccaH et al. 2000a). 

We conclude that an extremely steep IMF is not a viable 
alternative to explain the high values of the Mg indices in ellip- 
ticals. 

5.3. The effect of stellar evolutionary tracks 

As akeady stated, our standard SSP models are based on the 
stellar tracks by Cassisi et al. (1999). Since differences exist 
among different sets of tracks, it is interesting to check their 
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Fig. 5. Effect of stellar tracks on the Lick indices of 12 Gyr 
SSPs models with various total metallicities [Z/H]. Filled sym- 
bols: SSPs adopting the Cassisi tracks, open symbols: SSPs 
adopting the Padova tracks. 



impact on the model indices. To this aim we have computed 
SSP models with Worthey et al. FFs, but varying the input 
tracks. The Padova stellar tracks and isochrones as available 
on the Web have been used, specifically those by Fagotto et 
al. (1994) for metalHcities [Z/H] = -1.69 and -0.69, and 
those by Salasnich et al. (2000), with solar scaled abundance 
ratios, for metallicities [Z/H] = —0.4,0,0.35. The fuel con- 
sumption theorem is adopted also for these SSPs, following the 
method described in Maraston (1998; 2003). The results are 
shown in Fig. ^ which compares the Mgb and (Fe) indices 
of 12 Gyr SSP models as a function of the total metallicities 
[Z/H], as obtained with the Cassisi tracks (fiUed circles) and 
the Padova tracks (open circles). 

The use of the Padova tracks produces slightly stronger 
Mg6 indices, at metallicities [Z/H] between ^ —0.5 and so- 
lar. This effect is most likely due to the cooler temperatures 
of the Padova tracks along the Red Giant Branch with respect 
to the Cassisi tracks (see also Maraston 2003). A very modest 
impact is present also for the (Fe) index. Note that at metal- 
licities above solar the (Fe) obtained with the Padova tracks is 
stronger than that obtained with the Cassisi tracks, while the 
corresponding Mg6 indices are consistent. This implies that 
the discrepancy with ellipticals data is slightly larger when 
the Padova tracks are used. The differences are however very 
small, and the conclusions drawn from Fig. |l] are not affected 
by our use of a specific set of stellar tracks. 

Recently, isochrones and tracks with super-solar [a/Fe] ra- 
tios became available (Salasnich et al. 2000 with [a/Fc] = 
-1-0.3; see also Bergbusch & VandenBerg 2001; Kim et al. 
2002). These improved upon previous calculations that consi- 
dered only the effect on nuclear reaction rates (e.g. Salaris et 
al. 1993), while the updated tracks have also included the ef- 
fect of a-enhancement on the stellar opacities. Here we use the 
Salasnich's computations in order to check the impact of a- 
enhanced tracks on the final index values of SSP models. For 
consistency, we compare the indices based on the two Padova 
sets, with solar scaled and a-enhanced abundance ratio. This is 
shown in Fig. ^ for the illustrative case of 10 Gyr old SSPs|^. 
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We use here 10 Gyr to avoid interpolation on the Salasnich et al. 
isochrones 



Fig. 6. Effect of a-enhanced stellar tracks on Mg2 , Mg6, 
(Fe) and H/3 of 10 Gyr SSPs (Salpeter IMF) with various to- 
tal metallicities [Z/H]. Dotted lines: solar-scaled tracks. Solid 
lines: [a/Fe] = +0.3 tracks. Tracks are from Salasnich et 
al. (2000), the Worthey FFS are adopted. 



For the same total metallicity [Z/H], the Mg and Fe indices 
of SSPs based on a-enhanced stellar tracks (solid lines) are 
lower than those based on solar scaled stellar tracks (dotted 
lines). This happens because the a-enhanced tracks are bluer 
than the solar-scaled ones at fixed total metallicity (Salasnich et 
al. 2000), because of the lower stellar opacities. So, by increas- 
ing [a/Fe] at constant [Z/H] the Mg indices actually decrease! 

This rather counterintuitive behavior is a consequence of 
the fact that, at given total metallicity, the increase of the 
[a/Fc] ratio produces only a small increase of the a elements, 
and instead a large decrease of iron. In fact, in solar propor- 
tions the total metallicity is by far dominated by the a-elements 
(^ 74% by mass), the iron-peak elements amounting only to ^ 
8%, the residual being contributed by elements produced by the 
p, s, and r processes (Trager et al. 2000; TMB02).For example, 
compared to solar proportion in a mixture with [a/Fe]=H-0.3 
the a elements will contribute slightly more than ~ 74%, but 
the iron peak elements will be reduced by almost a factor of 
two, down to ^ 4%. Therefore, the main effect is to decrease 
iron rather than to increase, e.g., magnesium. Since iron is the 
most effective electron donor (e.g. Salasnich et al. 2000), the 
lower abundance of iron in enhanced [a/Fe] mixtures has the 
effect of decreasing their low-temperature opacities, which in 
turn determine an increase of the temperature of the RGB, and 
finally a decrease of the Mg indices because these are stronger 
in cool stars. 

The bluing of the isochrone also affects the H/3 line in- 
dices, which are stronger when a-enhanced tracks are adopted 
(Fig. ra, upper right panel). Therefore, in order to reproduce the 



C. Maraston et al.: Lick indices of Bulge stellar populations. II. Implications for SSP models and elliptical galaxies 



11 



observed Mgj, (Fe) and H/3 indices, the a-enhanced tracks 
require older ages. 

It should be noted that for the models of Fig. || the pa- 
rameter expressing the chemical composition in the Worthey 
et al. FFs (referred to as [Fe/H] in Worthey et al. 1994) has 
been considered to represent the total metallicity. As discussed 
in Sect. ||, this might be not entirely correct, since a variety of 
methods have been used, both photometric and spectroscopic, 
to determine such parameter for the stars of the Lick sample. 
If the [Fe/H]'s values are used, the metallic indices shown in 
Fig. ^ decrease even further. 

Fig. ^illustrates that accounting only for a-enhanced stel- 
lar tracks to compute a-enhanced SSP models, while using the 
same FFs, affects the indices in the wrong direction with re- 
spect to the locus occupied by the metal-rich GCs and ellipti- 
cals data in Fig. 0(Mg2 is affected more than (Fe)). A fully 
consistent exploration of the effect of the [a/Fe] ratio requires 
also the use of fitting functions depending on [a/Fe]. This is 
done in TMB02. 

5.4. Summary 

The conclusions of the previous paragraphs are the following. 

1 . Differences in the available sets of FFs do not affect the 
integrated Mg2 and (Fe) indices, due to the low contribution 
to the total continuum flux of those evolutionary phases where 
the FFs are mostly discrepant. 

2. For both Mg^ and (Fe) (and in general indices measured 
in the optical) the most important contributors are MS TO, 
RGB and HE stars. These are the evolutionary phases where 
the FFs need to be best constrained from stellar data. 

3. Results 1 and 2 hold for metallicities [Fe/H] ^ — 1, and 
are largely independent of the age and stellar tracks used. 

4. At very low metallicities (less than ZqIIO) the lower 
MS appears to dominate the value of the SSP Mg2 index. 
Uncertainties in the FFs and in the mass function jeopardize 
the calibration of the theoretical indices with the GC data. 

5. At metalUcities ^ ZqI2, the slope of the solar scaled 
Mg2 vs (Fe) relation for SSP models seems quite robust. One 
possibility to get strong Mg2 indices in combination with weak 
(Fe) without invoking a super-solar [a/Fe], is to adopt a very 
steep exponent for the IMF. However, IMF exponents as large 
as 4 5 have to be used in order to encompass the locus oc- 
cupied by ellipticals. Such values are ruled out by other con- 
straints. 

6. When a enhanced tracks are used the Mg2 and (Fe) in- 
dices become weaker, due to the blueing of the isochrone. 
Therefore, in order to match the observational data of GCs and 
ellipticals, abundance effects have to be accounted also in the 
FFs. 

6. Model Lick indices: comparison withi thie data 

In this section we compare quantitatively the indices for our 
sample of GCs with the models. Good models have to fulfill 
two requirements: i) the metallicities obtained using different 
line-strengths have to be consistent; ii) the metallicities derived 
from the models have to be in agreement with those determined 



independently from, e.g. spectroscopy of stars in GCs or CMD 
fitting. We already know from Fig. |l] that condition i) will not 
be fulfilled at [Z /H] ^ —0.6, because the GCs data deviate from 
the models. It is however important to check quantitatively the 
discrepancy. For a comparison of the standard SSP models used 
here with Magellanic Clouds clusters, we refer to Beasley et 
al. (2002). 

In this section we also explore the significance of other 
Lick indices as metallicity indicators. Finally we calibrate the 
B aimer lines. 

6. 1. Chemical compositions from Mg and Fe indices 

Fig. ^ compares the metallicities [Z/H] derived from the stan- 
dard SSP models, with those provided by the revised compi- 
lation of Harris (1996), which is largely based on the Zinn 
& West (1984) scale, for each GC of our sample. The model 
[Z/H] is obtained by interpolating the Mg2 index (left-hand 
panel) and the (Fe) index (central panel) on the SSP models 
(12 Gyr) based on the Cassisi tracks plus the FFs by: Worthey 
et al. (open circles), Buzzoni et al. (triangles), and Borges et 
al. (asterisks). The errorbars connect the minimum and maxi- 
mum model metallicities obtained by subtracting (adding) the 
observational errors to the measured values. For the empirical 
metallicities, a conservative error of H-0.2 dex has been con- 
sidered. The right-hand panel refers to a model (Fe) index 
as obtained by using as a metallicity input in the Worthey et 
al. FFs, the value of the iron abundance [Fe/H] and assuming 
an [a/Fe] = +0.3 (see Sect. Also in this case, the interpo- 
lation is with the model total metallicity [Z/H] . 
Fig. shows that for standard SSP models: 

i) the total metallicities [Z/H] as derived from the Mg2 in- 
dex (left-hand panel) are well consistent with the empirical 
scale of Zinn & West at metallicities lower than ~ —0.5. 
For NGC 6553 and NGC 6528 the models give [Z/H] values 
somewhat in excess of the values on the Zinn & West scale, 
but would agree with the near solar abundance indicated by 
the spectroscopic observations (Barbuy et al. 1999; Cohen et 
aL 1999). 

ii) the total metallicities [Z/H] as derived from the (Fe) in- 
dex, when the latter is computed with [Z/H] as input of the 
FFs (central panel), are systematically lower than the empirical 
ones, by roughly 0.3 dex. 

iii) the total metallicities [Z/H] as derived from the (Fe) in- 
dex, when the latter is computed by using instead the iron abun- 
dance [Fe/H] and assuming a 0.3 dex a-enhancement as input 
of the FFs (right-hand panel) is consistent with the Zinn & West 
values. 

On the basis of these evidences we conclude that: 
i) the standard SSP models, i.e. those based on the Milky 
Way calibrated FFs, do not underestimate the Mg2 index; rather 
they overestimate the (Fe) index. The disagreement between 
the models and the GC (and ellipticals) data would then point 
towards an iron deficiency, as opposed to magnesium enhance- 
ment, at virtually all metallicities. Suggestions in this direc- 
tion can be found in Buzzoni et al. (1994), Trager et al. (2000), 
TMB02; 
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Fig. 7. Comparison between the SSP-derived total metallicities from Mgj (left-hand panel) and (Fe) (central panel) for our 
sample of GCs, with the empirical metallicity scale [Fe/H] as compiled by Harris (1996), which is largely based on the Zinn & 
West (1984) scale. The standard SSPs, i.e. those based on MiUcy Way calibrated FFs, used here adopt the Cassisi tracks and: the 
Worthey et al. FFs (open circles), the Buzzoni et al. FFs (open triangles), the Borges et al. FFs (asterisks). The age of these SSPs 
is 12 Gyr. The SSP-derived [Z/H] is obtained by interpolating the observed indices on the model grid, separately. Diagonal lines 
show the 1 to 1 relations. In the right hand-panel the (Fe) index of the SSP models is derived by plugging [Fe/H] in the FFs 
instead of the total metallicity [Z/H](see Sect. 3), for the only case of the Worthey et al. FFs. 



ii) the model-derived total metallicities [Z/H] are in agree- 
ment with the metallicities on the Zinn & West scale which are 
refeiTed as to [Fe/H] (see Sect. 3.1). 

In the left-hand panel of Fig. ^ two of the four clusters at 

[Fe/H]zw 1.4 (NGC 6218 and NGC 6981) have a too 

low Mgj-derived metallicity, when the Worthey et al. or the 
Buzzoni et al. FFs are used. The SSPs with the Borges et al. 
FFs, instead reproduce the Mgj indices of these two specific 
clusters, because of the lower Mg2 indices of these FFs (see 
Sect. 4.2). Since at low Zs the Mg^ index is dominated by 
the lower MS component (Sect. 5.3), dynamical effects strip- 
ping low mass stars could be responsible for the observed low 
Mg2 indices of these particular objects. Piotto et al. (2001) 
show that indeed the mass function of NGC 6981 is consis- 
tent with a power-law with a slope flatter than the Salpeter one. 
We note that NGC 6218 is the object of our sample with the 
poorest sampled light (paper I). 

Fig. ^ compares the metallicities of each GCs as obtained 
from the Mgj (left-hand panel) and (Fe) (right-hand panel) in- 
dex using the TMB02 a-enhanced SSP models, with 12 Gyr 
and [a/Fe] — +0.3 (solid black lines in Fig. |l]). These models 
include the dependence of the fitting functions for /* on the 
[a/Fc] parameter. 

The total metallicities [Z/H] as derived from Mgj and 
(Fe) are in excellent agreement with each other. This results 
from having taken a super solar [a/Fe] abundance ratio in the 
models into account. The metallicities are in excellent agree- 
ment with those in the Zinn & West scale, over the whole 
range covered by our sample. The assumed [a/Fe] = +0.3 
for the GCs is in agreement with the [a/Fc] measured spec- 
troscopically (see e.g. Carney 1996; Salaris & Cassisi 1996). 
For the Bulge GC NGC 6528 recent spectroscopic abundance 
determinations give: Z ~ Z©, [a/Fe] - 0.3, [Fc/H] 0.3 



(Barbuy et al. 1999; Origlia et al. 2001). We acknowledge that 
some controversy exists in the literature: Carretta et al.(2001; 
see also Cohen et al. 1999) give [Fe/H] ~ +0.08 for 6528 and 

[Fe/H] 0.06 for NGC 6553, [a/Fe] - 0.2 + 0.4. It would 

be extremely important to pin down the element abundances for 
these objects which are the most metal-rich calibrators avail- 
able. 

CaiTetta & Gratton (1997) provide a metallicity scale for 
GCs which should reflect the [Fe/H] abundance. The right-hand 
panel of Fig. |^ compares the model-derived [Fe/H] with this 
scale. The metallicities in the Zinn & West scale are trans- 
formed into the CaiTetta & Gratton scale by adopting the re- 
lation provided by the authors. This comparison is rather poor. 
On this scale the iron abundance of our GC sample appears 
clustered on two values, while the observed (Fe) indices span 
a considerable range. It seems difficult to reconcile our index- 
based [Fe/H] with those in the Carretta & Gratton scale. 

We conclude this section discussing the case of 47 Tuc, 
which allows us to demonstrate the importance of comparing 
data and models which are set on the same system (see Sect. 
3.1, Table 1). If we use the indices measured in the Burstein et 
al. (1984) system, 47 Tuc has [Z/H] ~ -0.50; if we use the 
indices measured in the Worthey et al. one, we get [Z/H] ~ 
—0.78. This last value is in perfect agreement with the metal- 

Hcity of 47 Tuc in the Zinn & West scale ([Z/H] 0.76, from 

Harris 1996). In addition we notice that the metallicity derived 
for 47 Tuc using the indices in the Burstein et al. (1984) is as 
large as that of NGC 6356, while the comparison of the RGB 
ridge lines of these two clusters indicates that 47 Tuc is more 
metal-poor (Bica et al. 1994). For a critical discussion focused 
on 47 Tuc see Schiavon et al. (2002). 
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Fig. 8. The metallicities of the GCs as derived from the a-enhanced SSPs of TMB02, with age of 12 Gyr and [a/Fe] = +0.3. In 
the left-hand and central panels, the total metallicity [Z/H] as derived with the Mgj and the (Fe) index, respectively are compared 
with the values on the Zinn & West scale. In the right-hand panel the predicted [Fe/H] is compared to the empirical scale of iron 
abundances by Carretta & Gratton (1997). 

Table 3. The metallicities [Z/H] derived with Mgj and (Fe) of standard SSPs adopting the Worthey et al. FFs (col. 3-4), the 
Buzzoni et al. FFs (col. 5-6) and the Borges et al. FFs (see left-hand and central panel of Fig. The adopted SSP age is 12 Gyr. 
The last 2 lines report the calibration of 47 Tuc, using the original data by Covino et al. (1995) and those obtained by us on 
their spectra by adopting the W94 index definition (see 3.2, Table 1). In Col. 2 the empirical metallicities of the GCs as given by 
Harris (1996) are listed. 



name 


[Fe/H]zw 


[Z/H] 


WFFs 


[Z/H] 


SSP 
BFFs 


[Z/H]so^f,^, 






Mg2 


(Fe) 


Mg2 


(Fe) 


Mg2 


NGC 6981 


-1.40 


-1.91 


-1.83 


-1.87 


-1.73 


-1.46 


NGC 6637 


-0.71 


-0.73 


-1.05 


-0.76 


-0.94 


-0.54 


NGC 6356 


-0.50 


-0.51 


-0.91 


-0.55 


-0.81 


-0.35 


NGC 6284 


-1.32 


-1.27 


-1.58 


-1.29 


-1.41 


-1.07 


NGC 6626 


-1.45 


-1.31 


-1.48 


-1.33 


-1.31 


-1.12 


NGC 6441 


-0.53 


-0.53 


-0.80 


-0.57 


-0.70 


-0.36 


NGC 6218 


-1.48 


-1.80 


-1.71 


-1.77 


-1.57 


-1.35 


NGC 6624 


-0.56 


-0.45 


-0.88 


-0.59 


-0.78 


-0.39 


NGC 6388 


-0.60 


-0.81 


-0.84 


-0.84 


-0.74 


-0.61 


NGC 5927 


-0.37 


-0.17 


-0.65 


-0.23 


-0.56 


-0.07 


NGC 6553 


-0.34 


+0.05 


-0.27 


+0.01 


-0.29 


+0.11 


NGC 6528 


-0.17 


+0.09 


-0.35 


+0.06 


-0.35 


+0.14 


47 Tuc(B84) 


-0.76 


-0.50 


-0.66 


-0.53 


-0.58 


-0.33 


47 Tuc (W94) 


-0.76 


-0.78 


-1.07 


-0.81 


-0.96 


-0.59 



6.2. The other Lick indices as abundance indicators 



In this section we check the behaviour of the other Lick indices 
as abundance indicators, by analysing the correlations with 
Mg b and (Fe), which we have shown (Fig. 1) are likely to trace 
the a elements and the Fe abundances, respectively. To this 
purpose, we divide the indices into three groups: those which 
should be predominantly sensitive to a-elements (Mg]^,Mg2, 
etc.); those which should trace the Fe abundance (Fe4383; 
Fe5782; etc.); and the others (CNi; CN2; etc.), for which the 
case is less clear. The correlations are checked with respect to 
the model predictions. The evolutionary populations synthesis 
of Maraston (1998, see Sect. 3) has been updated for the com- 
putations of the whole set of Lick indices given in Worthey 
et al. (1994), plus the higher-order Balmer lines of Worthey 
& Ottaviani (1997). A detailed analysis of the contributions to 



these indices, similar to that given in Sect. 4 for Mgj and (Fe), 
will appear in Maraston (2003). 

Fig. H shows the correlation of the first group of indices 
with Mg b. The Mg indices: Mg b, Mg-^ and Mgj are very well 
consistent with each other and can be used as tracers of a- 
elements. The line-strength Ca4227 still appears to correlate 
with Mgb, although displaying slightly smaller values, while 
the other supposed Calcium-sensitive line Ca4455 does not. 
The reason for such mismatch is not clear to us. We suspect 
a calibration problem. The titanium-oxide indices TiOi and 
Ti02 behave consistently with Mg b at low metallicities, while 
at high metallicities these indices would overestimate the a- 
element abundance with respect to Mg b. However, TiO is con- 
tributed by M-type stars which may be poorly treated in the 
models and which are present only in the metal-rich clusters. 
The NaD index scales with Mg b although with some scatter. 
As discussed by TMB02, the reason for the scatter is most 
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Fig. 9. Calibration of a-sensitive Lick indices. Standard SSPs like in Fig. 1. Filled symbols denote our sample GCs, the large 
open symbol the average value of the Bulge light in the Baade window. 



probably the absorption by interstellar medium affecting ob- 
jects close to the galactic plane. Because of this possible source 
of contamination the NaD index is a problematic metallicity in- 
dicator for stellar populations (see also Burstein et al. 1984). 

Fig. |l^ shows the second group of indices that should trace 
the iron abundance. Indeed, all the plotted indices correlate 
very well with (Fe). The models fit well most of these indices. 
The only exception is Fe5782, which is underpredicted by the 
models, suggesting an offset in the FFs. This offset most prob- 
ably originates from the low signal-to-noise of the Lick/IDS 
data for Fe5782, which results from the extreme weakness of 
this line (S. Tmgev, private communication). 

It is worth noting that the C24668 index is tighly correlated 
to (Fe), although is supposed to trace the abundance of Carbon 
(see Trager et al. 1988). 

In Fig. [n] we plot the remaining indices vs Mg b. CNi and 
CN2 seem to trace the same elements as Mg b, though with an 
offset, especially at high metallicity and in the CNi. At least 
part of the effect could be due to the fitting functions not in- 
corporating stars belonging to high-metallicity GCs. As also 
discussed in Paper I, GCs have stronger CN indices compared 
to the bulge field (in Baade's Window) at the same value of 
(Fe). This may be caused by GCs stars accreting the CN-rich 
ejecta of AGB stars. The low-density environment of the field 



prevents a similar accretion on field stars. Therefore the models 
at high metallicity are lower than GC data because the fitting 
functions at high metallicity are contructed with field stars. 

Finally also the G-band seems to follow magnesium, but 
the relation is more scattered. 

In Fig.s H to |ll], the large open symbol shows the average 
index values of 15 Bulge fields, located in Baade Window. In 
all indices, the value of the Bulge field is consistent with an av- 
erage metallicity close to that of the most metal-rich GCs. This 
is quantitatively confirmed by the detailed metallicity distribu- 
tion of the bulge stars from optical-infrared color-magnitude 
diagrams (Zoccali et al. 2002). The large value of the NaD in- 
dex for the Bulge average field is again most probably due to 
contamination by interstellar medium, as discussed above. 

6.3. Balmer lines 

As well known Balmer lines are strongest in A-type stars, 
and become pregressively weaker for decreasing temperatures. 
In synthetic stellar populations their strength is sensitive to 
the temperature of the main sequence turnoff, hence to the 
age. Therefore the Balmer line strengths are used to estimate 
the age of e.g. elliptical galaxies, in an attempt at break- 
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Fig. 10. Model calibration. Iron-sensitive Lick indices. Model and data like in Fig. 
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Fig. 11. Model calibration of remaining Lick indices. Model and data like in Fig. |9[ 



ing the age/metallicity degeneracy (e.g. Worthey et al. 1992). 
However, turnoff stars are not the only potential contributor to 
the strength of the Balmer lines: Horizontal Branch (HB) stars 
may be as warm or even warmer than the turnoff. Actually, the 
H/3 index is perhaps more sensitive to the temperature distri- 



bution of the HB (the HB morphology) than to any other pa- 
rameter (Worthey 1992; Barbuy & de Freitas Pacheco 1995; 
Greggio 1997; Maraston & Thomas 2000). Therefore, attempt- 
ing to break the age/metallicity degeneracy with Balmer lines 
indices one runs into the age/HB morphology degeneracy. In 
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Fig. 12. Impact of the morphology of the Horizontal Branch on 
up. Models refer to SSP ages of 10 and 15 Gyr, and [Z/H] from 
—2.25 to 0.67 as labeled in the Fig.. Solid lines show models 
in which no mass loss has been accounted along the RGB. The 
dashed lines are 15 Gyr old models in which mass-loss has been 
applied at metallicities — 0.5, in order to reproduce the ob- 
served H(3 of metal poor GCs (see Maraston & Thomas 2000). 
The GCs of our sample are plotted according to their observed 
HB morphology, by means of the HBR parameter (from the 
Harris 1996 catalogue): pure blue HB (HBR=1, solid symbols); 
red HB (HBR=0, open symbols); intermediate morphologies 
(crosses). The open circle without errorbars shows the average 
value of 15 Bulge fields in Baade window. The very thick line at 
solar metallicity shows the locus of the solar metallicity mod- 
els, as corrected to take into account the bluening of the tracks 
due to a-enhancement (see Fig. 

the modeling the HB morphology cannot be derived from first 
principles, because of the role played by mass-loss on shaping 
the HB morphology . Additionally, possible dynamical effects 
(e.g. Fusi Pecci et al. 1993) may determine anomalous HB mor- 
phologies at a similar total cluster metallicity (the 2nd parame- 
ter problem in Milky Way GCs). Therefore the effect of the HB 
morphology to be used in SSP models needs to be calibrated on 
the Balmer indices. 

For our models this is done in Maraston & Thomas (2000), 
to which we refer for more details on the model parameter. In 
that work the mass-loss to be applied at every SSP metallicity 
was calibrated in order to reproduce the Hj3 line of the GC sam- 
ple by Burstein et al. (1984), Covino et al. (1995) and Trager et 
al. (1998). Here we check if those calibrated models are able to 
reproduce the Balmer lines measured for our sample sample. 

Fig. |l^ shows two sets of SSP models obtained with differ- 
ent prescriptions for the RGB mass loss, various metallicities 
and two ages (10 and 15 Gyr). Aimimg at calibrating the model 
Balmer fines with metallicity, we have chosen as x-axis the in- 



dex [MgFe] TMB02 show that this index by washing out 
[a/Fe] effects, is able to trace total metallicity. 

The solid lines connect models in which no mass loss is ap- 
plied to the RGB. In this case the morphology of the HB is red 
(i.e., all HB lifetime is spent on the red side of the RR-Lyrae 
location at logTcff ^ 3.85), except at the very low metallic- 
ity [Z/H] = -2.25, where the HB is spent at log Toff ^3.85 
even when no mass-loss is applied. The dashed line shows the 
15 Gyr SSPs with [Z/H] < — 0.5, in which mass-loss has been 
applied on the RGB according to canonical prescriptions. This 
leads to extended HBs, with blue morphologies (i.e. the whole 
HB lifetime is spent the blue side of the RR-Lyrae location) at 
[Z/H] = —2.25, and intermediate HB morphologies at metal- 
ficities between [Z/H] = -1.35 (~ 84% of the HB lifetime is 
spent blueward the RR-Lyrae) and [Z/H] = -0.55 (- 10% of 
HB lifetime is spent blueward the RR-Lyrae). 

The cluster data are plotted according to their observed HB 
morphologies, (B(lue)HB: filled symbols; R(ed)HB: open sym- 
bols; I(ntermediate)HB: asterisks) by means of the HBR pa- 
rameter (Harris 1996). 

The calibrated models by Maraston & Thomas (2000) are 
able to reproduce the observed Hf3 of our sample GCs. In 
particular they reproduce the relatively strong li/3 (~ 1.9 A) 
of NGC 6388 and NGC 6441, which are metal-rich clusters 
([Fc/H]zw = -0.6; -0.53) with an extension of the HB to the 
blue (Rich et al. 1997). The percentage of HB stars that is found 
blueward the RR Lyrae gap is: - 15% for NGC 6388 (Zoccafi 
et al. 2000b) and - 13% for NGC 6441 (M. Zoccali, private 
communication). This is well consistent with the HB evolution- 
ary timescales of the Maraston & Thomas (2000) models as 
given above. 

Concerning the most metal-rich objects NGC 6528 and 
NGC 6553, their relatively strong H/3 cannot be ascribed to 
HB effects since both clusters have a red Horizontal Branch 
(Ortolani et al. 1995). Part of the effect can be explained in 
terms of a-enhancement at high metallicities, without invoking 
young ages which would be in contradictions with CMD de- 
terminations (Ortolani et al. 1995). In Sect. 5 we have shown 
(Fig. ^ that a-enhanced tracks are bluer than the correspond- 
ing solar-scaled ones. As a consequence the H/3 lines are 
higher by 0.13 A at solar metallicity. The thick line in Fig. [T^ 
connects the SSP models as shifted by this amount. Note how- 
ever the rather large errorbar on the H/3 of NGC 6553. 

The index of the average light of the Bulge is shown as 
an open symbol without errorbars, and sits on the ~ 15 Gyr 
model with red HB. The low H/3 line does not leave room for 
intermediate age stars in our sampled Bulge fields. 

Finally we calibrate the higher-order Balmer lines 
HSa, a-fA, aSp and H7F (Worthey & Ottaviani (1997) of the 
same models, in Fig. |l^. It can be appreciated a very good con- 
sistency between the H/3and the higher-order Balmer lines. 

7. Summary and Conclusions 

Synthetic Lick indices (e.g., Mgj, (Fe), H/3, etc.) of stellar pop- 
ulation models (SSP) are calibrated over a range of metallici- 

' [MgFe]= v/Mgfc- (Fe> 
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Fig. 13. Calibration of the higher-order Balmer lines 
USa, USf, H7Aand H-fp (Worthey & Ottaviani (1997) 
of the standard SSP models, see Figure [T^. 



ties that extends up to solar metallicity using a sample of galac- 
tic globular clusters which includes high-metallicity clusters of 
the Galactic bulge. These data allow us to investigate empir- 
ically a well known property of elliptical galaxies, known as 
"magnesium overabundance" (Worthey et al. 1992), where the 
observed Mg indices in ellipticals are much stronger at given 
iron than predicted by standard models. By standard models 
we mean those constructed using stellar templates that are a- 
enhanced at low metallicity but assume solar elemental propor- 
tions at high metallicity. This effect has been generally inter- 
preted in terms of a-enhancement of elliptical galaxies, even 
if the bulk of their stellar population is metal rich. However, 
such conclusion rests on two assumptions: i) the models that 
are believed to represent solar scaled elemental ratios are cor- 
rect, and ii) the Lick Mg and Fe indices trace the abundance of 
the corresponding element. 

Using our GC database we have checked empirically both 
assumptions by comparing the Lick indices of the GCs spec- 
tra with those of ellipticals. The result is that the magnesium 
and iron indices of the metal-rich GCs, of the integrated light 
of the Galactic bulge, and of elliptical galaxies, define a fairly 
tight correlation in the Mg 6-(Fe)diagram, with elliptical galax- 
ies lying on the prolongation of the correlation established by 
the GCs, i.e., also the metal-rich GCs of the Milky Way bulge 
exhibit the "magnesium overabundance" syndrome. Since the 
GCs are indeed known to have enhanced [a/Fe] ratios from 
stellar spectroscopy ([a/Fe] ~ +0.3), we conclude that the 
interpretation of elliptical galaxy spectra in terms of "magne- 
sium overabundance" is indeed correct. The comparison with 
the GCs further allows us to point out that, rather than "magne- 
sium overabundance", the enhanced [a/Fe] ratio in GCs, and 
ellipticals is most probably the result of an iron deficiency with 



respect to the solar values (see Sect. 6.1). This agrees with pre- 
vious suggestions (Buzzoni et al. 1994; Trager et al. 2000b) and 
with the recent results by TMB02. The enhanced [a/Fc] ratio 
implies short formation timescales for the bulk stellar popula- 
tion, an important constraint for formation models of elliptical 
galaxies and bulges (e.g., Matteucci 1994; Thomas et al. 1999), 
which is difficult to reconcile with current semianalytic models 
of galaxy formation (Thomas & Kauffmann 1999). 

In parallel, the comparison of the SSP model with the GC 
data has allowed us to shed some light on the models them- 
selves. Around solar metallicity, the standard models are based 
on the stellar indices of Milky Way disk stars (by Worthey et 
al. 1994 or Buzzoni et al. 1992; 1994), and therefore would 
reproduce the indices of stellar populations with solar-scaled 
elemental proportions. This explains why they fail to repro- 
duce the Mg 6-(Fe)correlation followed by galactic GCs and 
the Galactic bulge, which are characterized by elemental ratios 
which are specific to the Galactic spheroid. This failure was 
also noted by Cohen et al. (1998) for both the metal rich globu- 
lars of the Milky Way as well as for those of M87, but was not 
attributed to an abundance effect. 

At metallicities [Z/H] < — 1 . the standard models use metal 
poor template stars that belong to the galactic halo, hence have 
supersolar [a/Fe] ratios. This explain why standard models 
successfully reproduce the Lick indices of the metal poor GCs, 
which also belong to the halo and are a-element enhanced. 
The conclusion is that the standard SSP models reflect abun- 
dance ratios which vary with metallicity. This clearly compli- 
cates their use as abundance indicators for extra-galactic steflar 
systems. 

We have then proceeded to compare the data to a new set 
of SSP models in which the [a/Fe] ratio is treated as an in- 
dependent variable (Thomas et al. 2002b). The result is that 
the Galactic GCs and bulge, as well as most ellipticals, are 
very well reproduced by coeval, old (12 Gyr) models with 
[a/Fe] — +0.3, and various metallicities. The uniqueness of 
this a-enhancement solution for the stellar populations of GCs, 
the bulge, and ellipticals was checked by thoroughly exploring 
the parameter space of the SSP models. We find that the Lick 
indices are little affected by the choice of the specific set of 
stellar evolutionary tracks or fitting functions. The only viable 
alternative to abundance effects, which can produce high values 
of the Mg indices coupled with low values of the Fe indices is 
a very steep IMF (much steeper than Salpeter). This solution, 
though formally acceptable, is practically ruled out by many 
other observational constraints for the clusters, as well as for 
the bulge and elliptical galaxies. 

A closer look to the Mg 6-(Fe)diagram reveals that ellipti- 
cal galaxies, unlike GCs, span a range of [a/Fe] values, from 
just marginally super-solar, to ~ +0.4. Since the Mg index cor- 
relates with the galaxy luminosity, the trend is in the direction 
of an increasing a-element overabundance with increasing lu- 
minosity (mass). Apparently, the more massive the galaxy, the 
shorter the duration of the star formation process (Thomas et 
al. 2002a). The origin of this trend remains to be understood. 

Our database of GCs was further used to check in an em- 
pirical fashion the effectiveness of the other Lick indices to 
trace the element abundances. Good indicators of a-elements 
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are found to be all the Mg lines (Mg2, Mg^ and Mg 6), 
and TiOi and Ti02 at subsolar metallicities. Also the index 
Ca4227 does correlate with the Mg indices, though a small off- 
set between the two might be present. Nearly all iron line in- 
dices (Fe4384, Fe4531, Fe5015, Fe5270, Fe5335) are found 
to display very tight relations against another. The indices 
CNi, CN2 and the G-band G4300 follow Mg. On the contrary, 
indices such as Ca4455, NaD and Fe5782 appear to be poorly 
calibrated, and we cannot recommend their use as abundance 
indicators for extra-galactic systems. 

The Balmer H/3 plus the higher-order lines by Worthey & 
Ottaviani (1997) H^^, HSp, H'^a, & H^f are very well repro- 
duced by the standard SSP models considered here (Maraston 
& Thomas, 2000; Maraston 2003), which indicates that they 
are only marginally affected by the [a /Fc] ratio (see Tripicco 
& Bell 1985; TMB02). Much more important for their correct 
modeling is to account for the Horizontal Branch morphol- 
ogy. In particular the rather high Balmer lines measured for 
NGC 6388 and NGC 6441 are modelled with a tail of warm 
Horizontal Branch stars ('^ 10% of the total HB population). 
These warm stars are indeed observed in the CMD of these 
two clusters (Rich et al. 1997), in a number {r-^ 10%, Zoccali 
et al. 2000b) which is in perfect agreement with the value re- 
quired to reproduce the strength of the Balmer lines. 

Finally, we point out that the Mg indices of very metal poor 
stellar populations ([Fc/H] ^ —1.8) are dominated by the con- 
tribution of the lower main sequence. Therefore, these indices 
are prone to be affected by the IMF and in GCs by the subse- 
quent evolution of the mass function due to the dynamical evo- 
lution of the clusters themselves. It follows that the Mg indices 
of very metal-poor stellar populations are not reliable metallic- 
ity indicators. 
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